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BRST symmetry has proven an invaluable tool in the perturbative quantisa- 
tion of gauge theories [1] so that its elevation to the non-perturbative level 
is clearly desirable. For one, BRST methods are crucial in the formulation 
of Schwinger-Dyson equations (SDEs) in covariant gauges, which seem gen- 
uinely non-perturbative, have been subject to study using various truncations 
at this level [2], and whose results are now subject to comparisons with cor- 
responding computations in Landau gauge from lattice gauge theory (see, for 
example, [3]). However, it has been demonstrated that a standard formulation 
of BRST symmetry invoked for lattice fields rigorously forces the partition 
function and path integrals of BRST invariant operators with a BRST invari- 
ant measure to vanish identically [4]. Rather than defining the configuration 
space in terms of some subset with no Gribov copies [5,6], one sums them 
all with alternating sign of the Faddeev-Popov determinant (only for small 
fields about = is this positive and thus the Jacobian of a change of vari- 
ables) and a complete cancellation takes place giving the nonsensical result 
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0/0 for lattice expectation values. This corresponds to attempting to resolve 
unity in the Faddeev-Popov trick via the partition function of a topological 
quantum field theory (TQFT) whose fields are the SU{n) group elements 
in the background of an external gauge or link field [7]. Such integrals of 
BRST invariant observables are topological invariants, the Euler character for 
SU{n) for the partition function, which all vanish giving the zeroes of the 
Neuberger problem. Equivariant gauge- fixing [8,9] evades this no-go theorem 
through a sequential gauge-fixing via coset space decomposition of SU{n), 
as in maximal Abelian gauge, so that the submanifolds of SU (n) implicit in 
this decomposition have non-zero Euler character. This formulation is how- 
ever distant from covariant gauge Schwinger-Dyson equations. Significantly 
though, there are quartic ghost couplings signalling that the damping term 
in the action for scalar Nakanishi-Lautrup auxiliary fields B is not BRST ex- 
act, B'^ 7^ s (something). This also signals the break down of the Neuberger 
argument. 

Quartic ghost couplings also arise in generalisations of the BRST and anti- 
BRST symmetry of Landau gauge through the so-called Curci-Ferrari (CF) 
"gauges" [10,11,12,13]. They also allow for a massive vector field while retain- 
ing BRST- (though not gauge-)invariance. But nilpotency of the (anti-)BRST 
algebra and thus unitarity are lost to be recovered, along with the original 
YM theory, in the massless limit. Nonetheless, this limit gives a theory which 
is local, covariant, BRST invariant, perturbatively renormalisable and close in 
spirit to Landau gauge. We will elucidate these details in light of the topo- 
logical approach to gauge-fixing and the Neuberger problem in the following, 
showing how these "bugs" possibly become features which may enable a non- 
perturbative definition of BRST. 

We assume antihermitian SU (n) generators T" so that ghost components C" 
are hermitian in order that C{x) — C"(x)T" = — On the lattice ghosts live 
on lattice sites Cj. In the continuum we deal with SU{n) gauge fields A^{x) ~ 
A°'J^x)T°' while on the lattice we denote the link fields as Uij for the link from 
lattice site i to j, [/,,,• = Pexp{J^^^ dz-A). The covariant derivative in the adjoint 
representation is D^- = ■ +[A^,-]. We also require a Nakanishi-Lautrup 
auxiliary 5'C/(n)-algebra valued hermitian field B{x) which lives on sites for 
the lattice theory, Bi. Note that the antihermiticity of the generators mean 
that the components of B are imaginary. Commutators and anti-commutators 
of variables will be separately indicated by [•, •] and {•,•}. 

For the continuum theory we have the BRST and anti-BRST algebras 



sAij_ — Df^C, sAfj_ — Di^C, 
sC^C^^ lie, C}, sC^C^^ 1{C, C} 



(1) 
(2) 
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sC^B-r\{C,C], sC^-B-r\{C,C]. (3) 

where for the lattice Eqs.(l) are replaced by 

sUij = CiUij — UijCj, sUij = CJJij — UijCj. (4) 

With these, the gauge-fixing is symmetric under C — > (7, C* — > —C. We have 
then sC + sC — {C, C} = which can be geometrically interpreted as the 
vanishing of a curvature in the extended space and so both ghosts and anti- 
ghosts arc Maurer-Cartan one-forms. We still have some freedom in specifying 
the variations of the auxiliary field B. We choose 



sB^m'C-^[B,C] + ^[C,{C,C}] 

sB^m'C -^[B,C]-^[C, {C, C}]. (5) 
An 5*^(2) group with generators cr"*", (J°, a~ specified by crM^ = a^B = and 



a-C = C, a^C = C, a+C = 0, a-C = 0, a°C = -C, a+C = C (6) 
enables us to replace nilpotency by the relations: 

2 2 + -2 2--I- 20 /rrX 

s = m a , s = —m a ,ss + ss — —m a . [7) 

The gauge-fixing of the Yang-Mills action can be achieved by the addition of 
a (anti-)BRST-invariant but not BRST-exact action, namely 

SGF^{ss-m^)W, (8) 



where for the continuum (respectively lattice) theory the most general choice 
for W is 

W^cont = ^Tr J d'xHA^y - CCC] , W,,t - ^ E Tr[Re[/,,- - ^QQ] (9) 

ij 



the first terms of which one recognises as the continuum (lattice) functional 
whose stationary points with respect to gauge transformations give the Landau 
gauge. 

The detailed form of the continuum Lagrangian after implementing the al- 
gebra has been given elsewhere [10,12]. We merely highlight specific terms. 
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Firstly, after integration by parts, s J d^xA^ = —2 / d'^xCd^ A^,. Acting with 
s then gives the ghost kinetic term and muhipher field term Bd^^A^, which 
are standard in covariant gauges. Secondly, the most complicated structures 
emerge from ss{CC). The damping term for the multiplier field emerges 
from here, also as in standard covariant gauges. But one ghost mass term, 
three-point couplings as well as quartic couplings (CC)^ and permutations 
thereof appear. Finally the m^{A? — ^CC) generates both a gluon and ghost 
mass term. 

Let us consider now some typical expectation value of a BRST invariant ob- 
servable of the link fields 0\U] in the lattice theory corresponding to the 
massive Curci-Ferrari gauge: 



where represents the auxihary fields, C, C and B. We shall now examine the 

precise relationship between this expectation value and that for lattice YM 
theory which, in terms of link variables, is well-defined even without gauge- 
fixing. We can factor into numerator and denominator the finite (on the lattice) 
integration over the gauge group, Vq = J Dg < oo and use the standard trick 
of exploiting the invariance of the measure, observable and YM action under 
gauge transformations to rewrite the expectation value as 



_ J DUe-'^^l^^O[U]Z[U] 
{0[U\}mCF - jjjUe-s.M[u\z[U] ^^^^ 



where 



Z[U] = j DgD(j)e-^°^^^'''t''^ (12) 



represents the partition function of a field theory in the group g and auxiliary 

variables (f) in the background of the link field U . But this is not a topologi- 
cal quantum field theory because the action of this theory is not based on a 
nilpotent algebra [14]. The consequence of this is that, unlike for TQFTs, this 
partition function depends on the background field, 7^ 0. The proof can 
be sketched as follows. Since the measure of Z\U] is independent of the hnk 
field the variation with respect to U acts directly onto the exponential of the 
action of the theory, Sgf\U^ ■,(!)]■, bringing the action into the measure. The 
variation 5/5U commutes with the operation ss — so that we effectively 
have 

/ DgD<t>{ss - m2)^e-^«-[^^'^l. (13) 



4 



But the integral of a BRST exact quantity with respect to an invariant measure 
still vanishes, despite the lack of nilpotency, thus 

^M=™^/fl,B^^^e-"P.,*, (14) 



which is not evidently constrained to vanish by any symmetry argument. Thus 
the partition function depends on the background link. Thus Z[U] cannot be 
factored out and cancelled except in the massless limit, so that YM expectation 
values can only be defined via the limit 

{0[U])ym^ ]im{0[U])rr^CF. (15) 
m— +0 



Unfortunately, the usual tricks of TQFT cannot help us in exphcitly evaluating 
Z[U] here: because of explicit dependence on U neither a semiclassical limit 
can be taken (as in TQFT) nor can the trivial link U = 1 he chosen (which 
for the lattice theory would result in some spin model [8]). We argue however 
that for generic m ^ the partition function Z[U] will be non- vanishing: 
the functional being introduced into the measure a la Faddeev-Popov trick 
is orbit- dependent and this is precisely what we require in order to lift the 
degeneracy between Gribov regions to avoid the Neuberger problem. We shall 
illustrate this below for the case of a simple one-link model. 

At any rate, we can give a final formula for the expectation value of a gauge- 
invariant observable in Yang-Mills theory in terms of the present construction: 



with Z[U] defined by Eq.(12). We can use the language of soft-meson theo- 
rems where the pion is a pseudo-Goldstone boson for massive quarks and thus 
describe Z[U] as the partition function of a pseudo-topological quantum field 
theory PTQFT. 

We can now elucidate how the Neuberger problem is avoided. Neuberger [4] 
considers the integral 

Io{t)= /D0e-^°-*^^O[t/,0]. (17) 



The measure D(f) is BRST invariant as is the action 5*0, which includes the 
Yang-Mills action. Expectation values in the theory are obtained for t = 1, 
namely (O) = /o(l)//i(l). However must also contain damping terms 
for the scalar B field integrations in Dcf) since even on the lattice these field 
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directions (unlike the link field U) are not compact. The damping term must 
be itself BRST invariant (upon an appropriate shift of the 5-field - actually in 
the standard case it is even BRST exact, = s (something) but this is not so 
relevant here) . Variation of Iq (t) with respect to t brings sF into the measure, 
the integrals of which vanish. Thus dio/dt — and Io{t) is t-independent. 
But for t = one has an integrand containing no ghost fields which vanishes 
by the rules of Grassmann integration: / DC = 0. Thus /o(0) = = /o(l) 
and all expectation values are of the form 0/0. This assumes I is well-defined 
at each step, which is only the case if Sq contains the damping term for B. For 
the massive-Curci-Ferrari case the structure of the action is different on two 
grounds: B"^ is not BRST exact, B'^ ^ s(something), but more importantly 
shifting B to B + ^{C,C} gives 

sTr = 2m^Tr(C6), (18) 



so that the damping term cannot be placed in 5*0 but must be placed in the 
term multiplied by t. There are two ways to do this, but keeping as close as 
possible to Neuberger's original argument we can reassign ssW tssW such 
that under derivation with respect to t a BRST exact term comes down into 
the measure. However, either way we cannot consider the t = limit as the 
functional integral then becomes undamped. This makes the Neuberger limit 
f — > fail, and so the usual proof fails. 

We now explicitly study this proposal in the context of the simple model 
introduced by Testa [15]. We consider an Abelian model with only two lattice 
sites, xi,X2 and thus only one degree of freedom, a link variable U which is 
parametrised through its phase U = e*"^. A is compact, A E [—^,^). A gauge 
transformation corresponds to shifting A by a difference {uj{x2) — uo{xi))/a 
which is a fixed quantity for any function uj[x). There are no plaquettes so the 
action is zero. The model essentially only contains topological information. 
The BRST and anti-BRST algebras for this simple gauge field theory can be 
written 



sA = C, sA = C 
sC = Q, sC = 

sC = i{B + CC), sC = i{~B + CC) 

sB = -i{m^ + B)C, sB = -i{m^ - B)C . (19) 
Note that B is now a real field, ghosts are Maurer-Cartan one-forms and 



s'C = m^C, = -m'C, (20) 

so nilpotency is lost. Using that sy[y4] = V'C and ^^[/l] = CV the gauge- 
fixing action for the massive Curci-Ferrari model here gives 
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S,f^{ss-m') {V[A]-^CC) 

= C[-V" + iV + 2^B + 2m^C]C + ^B^ + iBV' - mV (21) 

where V[A\ is constrained only by the requirement of periodicity under A — > 
A + 2Tx/a. Now this action appears Gaussian in all fields because quartic terms 
(7^(7^ all vanish since there is only one species of Grassmann field. Integration 
out of cither ghosts or scalar field B will upset this. In this case we can shift 
B via b — B + CC and avoid this. The action then becomes 

C[-V" + 2m^i\C + ihV + ih^ - mV (22) 

and integration out of h will give the Gaussian •v/7rexp(— ^(y[A])^)/-\/^. The 
procedure is dependent on the gauge parameter ^ for the same reasons as ^ is 
[/-dependent, but we will consider the case closest to the Landau gauge, ^ — 
for which we obtain then the delta function in the measure. Integrating 
out ghost fields gives for the partition function of the original theory: 

7r/a 

Z= 1 rfA(-y"[A])5[y'[A]]e™'^[^l (23) 

— 7r/a 

which should be compared with the partition function for the unfixed theory: 
ri^JA^^-Kja. 

We see that only stationary points of y\A\ contribute. They are weighted 
by the second derivative which would otherwise correspond to the Faddeev- 
Popov determinant; both signs of V" can appear. We crucially see the addi- 
tional orbit-dependent weighting exponential in the mass consistent with our 
observations above. If there are many such stationary points - Gribov copies 
- all of them will be summed over. But since y\A\ is not a gauge-invariant 
functional of A, they will generally come with different weight unless they 
represent degenerate stationary points of V\A\. For m = we recover the 
Neuberger pathology: critical points cancel according to the sign of the sec- 
ond derivative. The weight factor breaks this degeneracy so that the 
partition function will not vanish and the 0/0 problem disappears. 

Let us see this more explicitly. For simplicity we work in units of lattice spacing 
now (a = 1) and choose = | sin^ A so that V'\A\ = \ sin 2 A and V"\A\ = 
cos 2^4. The "gauge-fixed" configurations are thus 

^fixed = -7r,-|,0,| (24) 

for which V" has values 1, —1, 1, —1 respectively while the original V has values 
0, |, 0, ^. The partition function is trivially evaluated to be 
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Z = 2{e^ - 1) 



(25) 



Note that the constraint that observables O be BRST invariant in this case 
means that dO\A\jdA = 0, thus BRST invariant observables are just con- 
stants, O = c. All integrals of O are just multiples of the partition function, 
cZ. Thus the expectation value (0[74]) = c in the unfixed theory. Setting 
m = before taking the ratio reintroduces the pathological 0/0 result for 
the expectation value. Conversely, keeping m small but finite while taking the 
ratio for the expectation value allows factors to safely cancel in numerator 
and denominator, giving the correct result for the expectation value, namely 
the constant c. 

Whether such an elegant cancellation in expectation values of observables can 
take place in general, or subtleties of the m ^ limit need to be taken into 
consideration is an open question. For example, the vanishing of the partition 
function at m = can also be understood as the vanishing of the Witten index 
of the underlying TQFT with its supersymmetry. This leaves open the possi- 
bility of spontaneous breaking of BRST symmetry which would jeopardise the 
BRST cohomology construction of a physical state space. This aside, there is 
the question whether renormalisation effects can hinder the program. This is 
related to the other hmits which must follow a lattice calculation: the con- 
tinuum (a 0) and thermodynamic (L oo) limits. Note that for a finite 
lattice, the lattice gauge group is a simple product of SU{n) gauge groups 
per lattice site. Thus the original Neuberger zero corresponds to O'^"' °^ 
which goes over to a single zero in the a — > limit, related to the remaining 
global gauge symmetry of the torus at finite volume [8] . In order to avoid the 
Neuberger problem reappearing in our proposal it is essential that the Curci- 
Ferrari mass remain non-zero for finite L but vanishing a, particularly in light 
of renormalisation of the mass. For sufficiently large lattice volume, there will 
be no L dependence in the appropriate (and only UV sensitive) renormalisa- 
tion constants for m; for example [16,17] have computed these in continuum 
perturbation theory. Thus renormalisation will not introduce L-dependence 
into m. This suggests that it is safe to take m — either as fast as L — > oo 
or independently. But we stress that these considerations are only heuristic. 

We reiterate that unhke approaches seeking to isolate the Gribov or fundamen- 
tal modular region in the space of gauge fields, this approach takes all Gribov 
regions into account. There is the proposal that in the infinite volume limit 
configurations on the common boundary of Gribov and fundamental modular 
regions will dominate the ensemble averages of gauge invariant observables 
[18]. However such "dominant" fields will always be some subset of those 
configurations contributing at finite (L, a) thus representing no contradiction. 
This is also explicitly evident in the simple model. 

Open work includes careful examination of the m — > limit for the full theory. 
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We mention here that the violation of nilpotency of the BRST algebra results 
in negative norm states appearing in the physical Hilbert space [19,16] defined 
according to the Kugo-Ojima criterion [20]. This is one way of manifesting the 
loss of unitarity. These states do not belong to Kugo-Ojima quartets. In [19] 
one sees explicitly that in the massless limit they become zero norm ("daugh- 
ter") states. Thus from the point of view of the Hilbert space, the massless 
limit is smooth, as it is also for perturbative Green's functions. However, 
whether the same can be said for non-perturbative Green's functions is an 
open challenging question. There also remain technical challenges to imple- 
menting quartic ghost couplings in the lattice framework. The introduction of 
additional auxiliary bosonic fields, as in Nambu-Jona-Lasinio models, may be 
a way forward in this problem. 

To conclude, we have shown that the massive Curci-Ferrari model overcomes 
the Neuberger problem for elevating BRST symmetry to the non-perturbative 
level on a finite lattice. Beyond this first step, the verification that the mass- 
less, continuum, and thermodynamic limit of this procedure is the physically 
relevant theory faces a number of difficult challenges still. 
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